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Abstract 

Density functional theory (DFT) simulations corrected by the intramolecular Coulomb repulsion 
are performed to evaluate the vibrational inelastic electron tunneling spectroscopy (lETS) of O2 
molecules on Ag(llO). Semilocal DFT calculations predict a spinless adsorbed molecule, however 
the inclusion of the U leads to the polarization of the molecule by shifting a spin-polarized molecular 
orbital towards the Fermi level. This has an important implication in lETS, because a molecular 
resonance at the Fermi level can imply a decrease in conductance while in the off-resonance case, 
an increase in conductance is the expected lETS signal. We use the lowest-order expansion on 
the electron- vibration coupling, in order to evaluate the magnitude and spatial distribution of the 
inelastic signal. This allows us to reproduce the experimental data [Hahn et al., Phys. Rev. Lett. 
85, 1914 (2000)] in: (i) the negative conductance variation observed in the vibrational spectra 
of O2 along the [001] direction, (ii) the spatial distribution of the conductance changes recorded 
over the O2 molecule for the 0-0 stretch and the antisymmetric 02-Ag stretch vibrations, (iii) 
the absence of signal for the center-of-mass and hindered rotations modes, and (iv) the lack of 
lETS signal for the molecule chemisorbed along the [110] direction. Moreover, our results give 
us insight of the electronic and vibrational symmetries at play. The vibrational frequencies need 
to go beyond the harmonic approximation in order to be compared with the experimental ones, 
hence we present a Morse-potential fitting of the potential energy surface in order to evaluate 
accurate vibrational frequencies. The final lET spectra are evaluated with the help of the self- 
consistent Born Approximation and the effect of temperature and modulation- voltage broadening 
are explored. This ensemble of results reveals that the lETS of O2 is more complicated that a 
simple decrease of conductance and cannot be ascribed to the effect of a single orbital molecular 
resonance. 
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I. INTRODUCTION 



Single-molecule vibrational spectroscopy is now possible thanks to the development of 
the inelastic electron tunneling spectroscopy (lETS) operating with the scanning tunneling 
microscope (STM)^. The STM is thus conferred with the ability to chemically analyze sur- 
faces on the atomic scal^^. lETS has rapidly become a mature technique and multimode 
analysis on molecular substrates have been accomplishecP as well as analysis of extended vi- 
brations or phonons^. The characterization and control of single molecule reactivity through 
the excitation of a specific vibrational mode is a challenging area for the STM-based lETS 
techniqu^^. Within this context, the interaction of O2 with metals, which plays a central 
role in many technological processes, has served as a model system to explore and control 
the fundamentals of gas/surface reactivity. 

The chemisorption and dissociation of O2 on silver surfaces have been widely investi- 
gated in last decades in an attempt to understand the catalytic properties of silver, which 
are extensively exploited in industrial-oxidation processes. Molecular beam experiments 
provide detailed information on how reactivity depends on the collision energy, the molecu- 
lar rovibrational or electronic state, and the surface temperature and coverage. Combining 
this technique with either electron energy loss spectroscopy (EELS), infrared spectroscopy, 
or thermal desorption spectroscopy (TDS), it is also possible to extract information on 
the energetics ruling the elementary gas/surface processes -activation energies, atomic and 
molecular adsorption energies, reactive paths. Thanks to these kind of studies we know 
that O2 dissociation on fiat Ag surfaces is characterized by rather large activation ener- 
gie^^ and, as a consequence only molecular adsorption is possible at crystal temperatures 
below 150 Kf^SHm^ ^j^e chemisorption of O2 on the Ag(llO) surface has been particularly 
controversial because of the initial disagreement regarding the nature and orientation of the 
chemisorbed molecule. Density functional calculations performed by Gravil et al^ shed light 
to this controversy showing the existence of two distinct chemisorption states with essen- 
tially equal energies. This theoretical finding was subsequently corroborated by EELS and 
TDS studies^^, though the final confirmation came with the STM-IETS investigations by 
Hahn et al]^^^. Compared with other spectroscopies, the great advantage of lETS operated 
with STM is that it provides simultaneous topographical and spectroscopical images with 
atomic resolution. This permits an almost direct identification of the molecular state. 
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Molecular O2 chemisorbs on Ag(llO) parallel to the surface over the hollow site on 
two possible configurations, one with the molecular axis along the [001] surface direction 
(O2[001]) and the other one along the [lIO] direction (O2[li0]). Single molecule vibrations 
are detected with STM-IETS on the O2[001] only. The recorded inelastic signal is a decrease 
in conductance for one of the modes (the 0-0 stretch) and it can be a decrease or an in- 
crease in conductance depending on the exact position of the STM tip over the molecule 
for the other detected mode (the antisymmetric 02-Ag stretch)^. The reasons behind this 
rich and complex lETS structure have been a matter of controversy that have recently been 
solvecP^. 

In lETS, the changes of conductance are recorded as a function of the tip-substrate volt- 
age. It is generally admitted that in tunneling, the excitation of a vibration mode leads to 
opening an inelastic channel for conduction, hence increasing the conduct anc#^^. However, 
early model calculations predicted that when a molecular resonance overlaps the Fermi level 
of the substrate, lETS could also give rise to decreases in conductance^^. Up to now, the 
02(001] is the only example, where such decreases in conductance have been observed exper- 
imentally. Strikingly, state-of-the-art lETS simulations, which use the electronic structure 
calculated with DFT as the initial ingredient, have been unable to reproduce the experi- 
mental conductance decreases. The simulations give conductance increases because in the 
DFT electronic structure the Tig orbitals are not at resonance with the Fermi leve F^ ' ^^ l This 
discrepancy led to speculate that the theoretical Fermi level was wrongly positioned. How- 
ever, questioning the position of the Fermi level and, thus, the molecular charge state in this 
system is at variance with the excellent characterization of the two experimental chemisorp- 
tion states given by this theorj^. Moreover, the Fermi level fitting, thought to force the 
TTg resonance, cannot reproduce the rich lETS structure of the experimental data showing 
increases and decreases according to the vibrational mode and the tip localization. 

In a recent letteil^, it has been shown that correcting the DFT electronic structure to 
adequately incorporate the on-site Coulomb repulsion (the so called DFT+U approach), the 
experimental conductance changes recorded over the 0-0 stretch mode and the antisymmet- 
ric 02^Ag mode are correctly reproduced. Furthermore, the DFT+U calculations permit us 
to gain extra insight on the actual electronic structure of this molecular system. Contrary 
to the previous speculations, rather than shifting towards the Fermi level, the TTg resonance 
follows a St oner-like splitting of its spin components, leading to the partial filling of one of 
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the spin components of the tt^ orbitals and, thus, to the spin polarization of the isolated 
O2 molecule on Ag (110). In the present contribution we present a comparative analysis of 
the DFT+U electronic structure of the O2[001] and ©2(1 10] states, discussing the properties 
that cause the respective survival and quenching of the intrinsic O2 paramagnetism. Our 
previous simulations of the conductance changes for the O2[001] state are extended to all 
the relevant molecule-surface modes and also to the ©2(1 10] molecule. In particular, the 
reasons for the existence or absence of lETS signal are explained in terms of the electronic 
and vibrational mode symmetries. There are two symmetry planes in the O2/Ag(110) sys- 
tem. We show that the possibility to excite a vibrational mode is determined by the overall 
symmetry character of the coupling states and mode respect to those two planes. However, 
the value of the lETS signal depends ultimately on the energies of the orbitals and on the 
strength of the electron vibration coupling. Also, we evaluate the lETS spectra with the 
help of the self-consistent Born Approximation and explore the effects due to temperature 
and modulation-voltage broadening. 

The reminder of the paper is organized as follows. Section |II A| describes the details of 
the DFT+[/ calculations performed to obtain the electronic spectra of O2 chemisorbed on 



Ag(llO). In Sec. II B, we discuss the necessity to go beyond the harmonic approximation 
in order to obtain accurate vibrational frequencies of the molecule-surface modes. The 



theoretical model used for the lETS simulations is explained in Sec. |IIC[ The calculated 
electronic structure is shown in Sec. [IITl The corresponding STM-IETS simulations are 
discussed in Sec.jTVj In Sec.[V[ we apply the SCBA model to calculate the lETS. Conclusions 
and final remarks are presented in Sec. |VI[ 



II. METHODS 



A. Details of the BFT+U calculations 

Density functional theory calculations are performed with the VASP codd^^l using plane 
waves with an energy cutoff of 515 eV and the projector augmented wave methocP^. We use 
a periodic supercell consisting of a six-layer slab separated by 10.45 A of vacuum and a (3x4) 
surface unit cell large enough to significantly reduce interactions among oxygen molecules in 
neighboring cells. The surface Brillouin zone is sampled with a 4x4x1 Monkhorst-Pack grid 
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of special k points. The exchange and correlation energy is calculated within generalized 
gradient approximation (GGA) using the PW91 functionaP^. Corrections to the on-site 
Coulomb repulsion are applied to the oxygen 2p electrons following the rotationally invariant 
DFT+U scheme proposed by Dudarev et al^^ as implemented in VASP. The screened on-site 
Coulomb interaction U is calculated from first principles as the energy cost for adding extra 
charge to O2 when adsorbed on Ag(llO). The value obtained with the method described in 
Ref. is close to 4 eV. However, we also verify that the results shown in this work, are 
reproduced for U varying within the range 2-4 eV. Such a verification is necessary because 
the calculated U value actually depends on the projectors and the used constrained method 
to fix the final electronic structure. The adsorption positions are optimized by fully relaxing 
the O atoms and the two uppermost silver layers until atomic forces are less than 0.02 eV/A. 

B. Vibrational frequency evaluation: beyond the harmonic approximation 

The diagonalization of the dynamical matrix gives the vibrational modes within the har- 
monic approximation. The method implemented in VASP uses finite differences to compute 
the second derivative of the hamiltonian with respect to a constant displacement AR of 
every coordinate of the active atoms, here the oxygen ones. In Table [T] we give the first 4 
molecular modes on the surface for atomic displacements, AR = 0.025, 0.03 and 0.04 A. 
The displacement of 0.04 A is the one that has the better overall agreement with the ex- 
perimental measurements. The dispersion is still considerable, despite DFT is known to 
yield frequencies within ^ 5% the experimental values^^. The use of a rigid displacement to 
estimate the dynamical matrix is too simple to capture the complexity of molecular modes 
on surfaces. Indeed, a small displacement will be probing a region too fiat of the potential 
energy surface (PES), and a large value will be probably in the anharmonic part of the 
PES. The situation becomes more difficult when we consider modes that span several tens 
of meV as is the present case. A AR that probes the harmonic potential corresponding to 
the center-of-mass coordinates (at ^ 30 meV) will be too large to probe the 0-0 stretch 
mode (at ~ 80 meV). 

The vibrational frequencies are accurately calculated if several points of the PES are 
obtained and fitted to a Morse potential. This is the strategy followed here. We use the 
eigenvectors of the simple finite difference method. These eigenvectors change little from 
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displacement to displacement (less than the frequencies do) and are a good starting point to 
plot the PES along a given mode. Figure [T] shows the PES along the 0-0 stretch mode and 
the corresponding Morse-potential fit. Solving the one dimensional the Schrodinger equation 
for this PES, we find a frequency of 83.626 meV. This value is now in good agreement with 
the experimental frequency of 82.0 meV. Therefore, we safely conclude that the uncertainty 
in the frequency evaluation is due to the anharmonicities of the PES and not to the DFT 
calculation. 

C. lETS simulations 

The lETS simulations of O2 on Ag(llO) are performed using two different methods. 
For the quantitative evaluation of lETS and its spatial distribution we use the many-body 
extension of the Tersoff-Hamman theory for the STlVp^^l. The current implementation is 
based on DFT results regarding the electronic structure and the electron- vibration couplings. 
Briefiy, the conductance change due to the excitation of a localized vibration to the lowest 
order in the electron-vibration coupling, v is given by the opening of the inelastic channel, 
leading to an increase of conductance, A(7^^e(r), and the change in the elastic channel to the 
same order in A(7e/a(r), which means a decrease in conductance. In Ref. [31j, an all-order 
theory shows that the decrease of conductance is related to the appearance of vibrational 
side-bands in the elastic electron transmission. The decrease is maximum when a molecular 
level is resonant with the Fermi energy and such that side bands appear at -^hVt. 

Due to the strong molecular character at the Fermi energy, approximations replacing 
the energy dependence by the behavior at the Fermi energy are not correct any longer. 
Hence, we use the full electronic structure and energy dependence for the finite bias of the 
measurements. In this case, we extend the customary lETS treatment! ^^ ' ^^ ' to new equations 
that keep all the energetic dependence. Hence, the inelastic contribution to the conductance 
change is given by: 

define _ 1 

a p(ro,^F + ey) 

2 

(l-/Kk))5(^F + ey-/if^-e,,k), (1) 
where p(ro, Ep + eV) is the local density of states evaluated at the STM tip's center Fq and 
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at the Fermi energy plus the energy of the electron at the corresponding bias, V. The Fermi 
distribution function is given by /(e) and Q is the frequency of the considered mode. Since 
we are using periodic-boundary calculations, the electronic states are Bloch states with band 
indexes n and m and k-vector k. Note that the local electron-vibration potential v couples 
states with the same k-vector. The elastic term is 



p{ro,EF + eV) 



E 



J]] (1 - /(em,k)) ('0m,kb|'0n,k) ^m,k(ro) S{em,k - hQ - e^,k) 
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(2) 



The factor (1 — f{en^u))S{EF + eV — HQ — e^^k) is responsible for the temperature-dependent 
onset of the vibrational signal. Indeed, this factor is strictly zero at zero temperature, if 
eV < hbarQ. In the evaluation of these equations the k-point sampling is critical because a 
large number of electronic states is needed to ensure numerical accuracy. The 5-functions are 
approximated by Gaussian functions of broadening smaller than the used bias. Typically, 
we have performed the simulations dit V = 150 mV with a gaussian broadening of 50 meV, 
except for the 0^0 stretch mode that is calculated at 1/ = 200 meV. The experimental 
change in conductance is simulated as the sum of the inelastic and elastic contributions. 

The second method is the use of the self-consistent Born approximation (SCBA)P2i3^ 
based on a parameterized hamiltonian following Ref. [31j. In spite of the hamiltonian sim- 
plification, SCBA is numerically difficult and the evaluation of the spatial distribution of the 
SCBA signal is beyond our capabilities, instead we have fixed the STM-substrate symmetry 
and evaluated the lETS as a function of the applied bias. Results are shown depending on 
temperature and modulation voltage of the experimental lock-in amplifier in Sec. [V| 



III. ELECTRONIC STRUCTURE OF O2 ON Ag(llO) 

Despite local and semilocal DFT have proven to be accurate in determining variational 
ground state properties of the many-body system such as the electron density, it cannot 
assure a reliable treatment of the spectroscopic properties. In many cases, the electronic 
structure is incorrect due to the inadequate description of the on-site Coulomb interaction 
that favors fractional occupancies of energetically close states. This is precisely what happens 
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with the electronic structure of the O2[001] chemisorption state calculated with GGA. In 
Fig. [2] (a) the projection of the density of states (PDOS) on to the O2 molecular orbitals 
(MOs) reflects the partial occupation of the four tt^ orbitals, despite the adsorbed molecule 
is capturing less than two electrons. These unrealistic populations are possibly at the origin 
of the discrepancy between the experimental lETS and the GGA-based simulations. To 
solve this failure we have performed DFT+U calculations of O2 chemisorbed on Ag(llO). 
The DFT+U is aimed to correct the electronic structure without perturbing the physical 
magnitudes for which the GGA excels. A detail comparison of these quantities is shown 
in Table [Til for the two chemisorption states. The good agreement obtained with the GGA 
values ensures the correct description of the system. 

The more realistic DFT+U electronic structure shows that the paramagnetic nature of O2 
is partially preserved when chemisorbed along the [001] direction. As shown in Fig.[2](b), the 
adequate description of the intramolecular Coulomb interaction breaks the spin degeneracy 
imposed by GGA and allows the full occupation of the two parallel-to-the-surface TTg orbitals 
(ttI) and the spin-up perpendicular-to-the-surface Tig (vrj-). Still the spin-down irj- is partially 
occupied in order to preserve the correct electron density ascribed to the molecule. As a 
result, this molecular orbital is at resonance with the Fermi level. The electronic structure 
of this chemisorption state can be viewed as a Stoner-like process in which the molecule 
reduces its total energy by splitting the orbitals according to their spin polarization. One 
becomes completely filled at the expense of the other stabilizing the full system. 

In contrast, no spin-splitting of the vrj- orbital is obtained for the O2[110] molecule. At 
first sight there are no significant differences on the structural properties (Table [ll]) and the 
GGA electronic structures of the two chemisorption states [compare Figs. [2] (a) and (c)]. 
These similarities are indeed consistent with the experimental observations of almost iden- 
tical chemisorption energie^ and equal probability to find any of the two configurations^. 
However, the small differences found in the GGA calculations already point to a slightly 
larger charge transfer from the surface to the O2[li0] molecule, which is also closer to the 
surface and more stretchecP^. Similar results were also found in Ref. [8j using a 3 x 2 sur- 
face cell. Experimentally, this picture agrees with the existence of two energetically-close 
modes that are attributed to the 0-0 stretching vibration of the O2[li0] (79.5 meV) and 
the O2[001] (85 meV)'^. Our DFT+U calculations show that the induced density for O2[110] 
is larger and the Coulomb repulsion less important. Interestingly, what the DFT+U results 
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highlight is a new difference between the two chemisorption states, namely, the paramagnetic 
character that is only preserved in the O2[001] molecule. 

The DFT+U electronic structures are used to generate the constant current simulations 
of Fig. [sj The experimental topographical images of the OafOOlpiand the OafllOP are well 
reproduced in both cases. The main difference between both states are the two protrusions 
appearing along the molecular axis in the former but not in the latter. As was demonstrated 
by the non-spin polarized DFT calculations of Olsson et al^^ ^ these protrusions are the 
fingerprints of the vrj- orbital. The absence of these features in the LDOS image of the 
O2[110] molecule is also understood on the basis of the stronger screening exerted on this 
molecule by the surface electrons. Note, finally, that the similar STM simulations obtained 
with non-spin polarized DFT and spin polarized DFT+U remark that information about the 
spin moment cannot be directly extracted from the STM topography. In next section, we 
will see that STM-based lETS can alternatively be used to identify this peculiar electronic 
structure in which a spin-polarized molecular orbital (MO) is at resonance with the Fermi 
level. 

IV. SPATIAL DISTRIBUTION OF CONSTANT CURRENT lETS SIGNALS 

Our STM-IETS simulations, based on the DFT+U electronic structure, are able to re- 
produce qualitatively, and many times quantitatively, the reported experimental findings: 
(i) the negative conductance changes observed in the vibrational spectra of the O2[001] 
molecule, (ii) the spatial distribution of the conductance changes recorded over the O2[001] 
molecule for the 0-0 stretch and the antisymmetric 02-Ag stretch vibrations, (iii) the ab- 
sence of signal for the center-of-mass and hindered rotation modes, and (iv) the lack of lETS 
signal for the molecule chemisorbed along the [110] directioiP^l. 

We start by analyzing the spatial distribution of the conductance change for the O2[001] 
molecule. A schematic diagram of the vibrational modes is depicted in Fig. [5| Figures [i] (a) 
and (b) show that the symmetrical 0-0 stretch mode leads to a decrease of conductance 
over the adsorbed molecule. In agreement with the experiment, the maximum absolute value 
of conductance change is displaced from the center of the oxygen atoms in the same way as 
the constant current STM image. This is reminiscent of the nj- orbital, showing that the 
conductance decrease due to this vibrational mode can be understood within the framework 
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of a single level resonant with the Fermi energy. Model calculations^^ demonstrated that the 
presence of a substantial density of states at the Fermi energy of large molecular character 
leads to an important lETS signal of negative sign, i. e. a reduction of conductance over 
the vibrational threshold, as observed here. 

The antisymmetric stretch mode is more interesting. Here, the oxygen atoms move 
vertical to the surface with opposite phases. Again, in good agreement with the experiment, 
we find that while the change in conductance is negative away from the molecular center, 
it is positive at the molecular center [Figs, [i] (c) and (d)]. In contrast to the 0-0 stretch 
mode, this behavior implies that the lETS signal involve the coupling with various molecular 
orbitals and, therefore, we should go beyond the single-level model to understand this lETS 
image. 

Simulations performed for the center-of-mass mode show a negative conductance change 
image rather similar to the one obtained for the 0-0 stretch mode. The maximum con- 
ductance decrease however is a factor of 5 lower in the center-of-mass, what explains the 
absence of signal in the experiments for this mode. The simulated efficiency for the hindered 
rotation mode, in which the atoms try to move opposite to each other along the [001] sur- 
face direction, is even smaller: less than 0.8% of conductance change. Neither is this mode 
observed in the recorded lETS data. 

All the above results can be understood with the help of Eqs. ^ and Q. As the nj- is 
the only MO at resonance with the Fermi level, the topography of the conductance changes 
is determined by the molecular orbitals coupling with such single resonance, i.e., by those 
MOs contributing with nonzero matrix elements in Eqs. ([T]) and ([2]). The symmetry analysis 
summarized in the table of Fig. [5] provides a preliminary idea of the MOs that can couple 
with the nj- for each vibrational mode. The weight of these couplings in the lETS image 
will ultimately depend on the orbitals energies and on the strength of the electron- vibration 
coupling as discussed below. 

There are two well defined symmetry planes in this system: the plane perpendicular to 
the molecular axis ay and the plane normal to the surface that contains the molecular axis 
a^. The regions where the wave functions take positive and negative values are schematically 
depicted for each MO in the upper panel of Fig. [5] with blue and red colors, respectively. As 
written in the table, the vrj- is antisymmetric (A) respect to and symmetric (S) respect to 
a^. Since the 0-0 stretch mode is symmetric respect to both planes, the coupling electronic 
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states should respectively be A and S. These conditions are only satisfied by the vrj- and 
the au orbitals. The same conclusion is obtained for the center-of-mass mode, which is also 
symmetric respect to both planes. The A and S characters of the antisymmetric 02-Ag 
stretch mode respect to and a^, force the coupling of nj- with the ag and the tt^ orbitals. 
Finally, the excitation of the hindered rotation mode, which presents a fully antisymmetric 
character, couples the nj- with the ttI orbital only. 

Once we know which orbitals are coupled by each vibrational mode, the topography of 
the lETS signal is better understood if the conductance change is divided in its inelastic 
contribution due to the opening of a new channel, hence positive, and the renormalization of 
the elastic channel due to the presence of the vibration that leads to a negative contribution. 
Figures [g] and [T] show the inelastic (left panels) and elastic (right panels) contributions for 
the four vibrational modes. According to Eqs. ([T]) and ([2]), the main diflFerence between 
them despite their sign, is that the inelastic part involves all possible intermediate electronic 
states, while the states contributing to the elastic part are restricted by the extra 5-function 
to an energy interval of the order of the mode frequency HQ about the Fermi energy. This 
strict energy restriction explains that in all cases the vrj-, being the only orbital at resonance 
with the Fermi level, dominates the topography of the elastic contribution. 

In contrast, the inelastic image can be composed by all the coupling states that according 
to the symmetry analysis contribute with nonzero matrix elements in Eq. ([T]) . The weight of 
each coupling state in the lETS image is partially determined by the denominator of Eq. ([T]) , 
penalizing the orbitals that are energetically distant from E^. The effect of the diflFerent 
weighting is refiected in the simulations of the inelastic contributions. Figs. [6] (a) and[7| (a) 
show that the inelastic contribution is due to the (at resonance with Ep) and not to 
the energetically distant in the 0-0 stretch and center-of-mass modes. Otherwise, there 
would be a certain intensity modulation along the molecular axis caused by the nodes of 
the au orbital. Similarly, the absence of nodes along the molecular axis and the maximum 
intensity obtained at the center of the molecule in Fig. [g] (c) show that the image is basically 
due to the tt^ and not to the ag in the antisymmetric stretch mode. In the hindered 
rotation mode. Fig. [t] (c) shows that the inelastic contribution is uniquely formed by the ni 
as predicted by the previous symmetry analysis. 

Now we have all the ingredients necessary to understand the lETS data recorded for 
the O2[001] molecule. In the 0-0 stretch mode, the elastic and inelastic contributions are 
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only due to the Trj- orbital. As the elastic part is much larger than the inelastic one, the 
excitation of this vibration leads to an overall reduction in conduction except between the 
O atoms where the tt^ orbital presents a nodal plane. Compared with the 0-0 stretch 
mode, the smaller elastic and inelastic contributions obtained in the center-of-mass mode 
are due to the weaker electron-phonon coupling, as can be roughly inferred from the different 
vibrational energy of each mode (see table |l]). In consequence, no lETS signal is observed in 
the center-of-mass mode, despite its symmetry similarities with the 0-0 stretch mode. The 
antisymmetric 02-Ag mode is more complex because involve different orbitals: the inelastic 
part is dominated by tt^ and the elastic part by the vrj-. The rich structure observed in 
this mode is a consequence of the different nodal structure of the two orbitals. The total 
conductance change is positive in the region between the two oxygen atoms because the 
elastic contribution is zero whereas the inelastic one is maximal. On the contrary, the 
conductance change is negative out of the internuclear region because there the elastic part 
is maximal and the inelastic part minimal. Finally, the absence of lETS signal for the 
hindered rotation mode is a consequence of two factors: the low DOS with vrli character 
around Ep and the weak electron-phonon coupling of this mode. 

We also perform simulations of the conductance changes for the O2[110] chemisorption 
state. The absence of experimental lETS signal is in agreement with our calculations that 
show changes of about 1% at most. The reason behind this low conductance changes is 
clear: as the vrj- is practically occupied, the probability to excite any vibrational mode by 
tunneling trough this MO is small. 

V. BIAS DEPENDENCE OF lETS: A SELF-CONSISTENT BORN APPROXIMA- 
TION STUDY 

In order to compute the bias dependence of the inelastic signal, we use the self-consistent 
Born approximation SCBA. Briefly, the SCBA consists in using the Born expression for 
the electron- vibration self-energy (lowest order in the electron- vibration matrix element) re- 
placing the free-electron Green's function by the full Green's functioii^. This last function 
is unknown, hence an iterative process is used where the starting step is the free-electron 
Green's function, and the second step uses the self-energy computed with the Green's func- 
tion obtained from the first step. This self-consistent procedure is stable and has the virtue 



13 



of generating high-order terms in the electron-vibration couphng. 

The SCBA has been implemented in DFT based approached^. But it is a costly method 
for evaluating STM scans as the ones shown herein. It is however feasible for a single STM 
conformation, where the bias is ramped and a (fl /dV^ spectrum is obtained. Here, we have 
simplified the problem and taken the tight-binding parametrization of Ref. ,3Jj 

As discussed in Ref. |43l the SCBA fails in the limit of large electron- vibration coupling. 
The main reason is because the SCBA fails in taking into account the correct contribution 
of the number of excited modes. Indeed, H. Nesd^ shows that SCBA yields the same result 
as an exact calculation if the electron-vibration matrix elements did not depend on the 
number of excited modes. If the average number of excited modes is in the order of one, 
then the SCBA gives a very good account of the vibrational excitation even in the case 
of extraordinary vibronic eflFects^^. Approaches going beyond SCBA are availabld^™ and 
are necessary for exploring realistic situation where the electron- vibration coupling is strong 
enough to excite several overtones. 

Figure [S] shows the results of the SCBA for the 0-0 stretch mode. Panel (a) shows the 
conductance, mainly given by the molecular resonance about the Fermi energy, and two 
sharp Stokes and anti-Stokes changes in conductance. In (b) the conductance derivative or 
(fl/dV^ is shown. In both cases three different temperatures are used, 6 K, 14 K and 20 
K. At low temperatures the d'^I/dV'^ shows a Fano profile as has already been discussed for 
the excitation of vibrations by electronic currentd^, however as the temperature increases, 
the peak becomes broadened and closer to a simple Lorentzian peak. It is the modulation 
voltage used in the lock-in detection that finally broadens and removes any particular detail 
from the excitation peak. Figure [s] (b) shows in full line the peaks after convolution^^ using 
a modulation voltage of 7 mV rms as in the experiment!^. 

VI. CONCLUSIONS 

We have conducted STM-IETS simulations, based on the DFT+[/ electronic structure 
and applying an extension of the Tersoff-Hamman theory to the presence of inelastic ef- 
fect^^^. Our results are able to reproduce the experimental lETS recorded by Hahn et 
al^^^. (i) the negative conductance changes observed in the vibrational spectra of the 
O2[001] molecule, (ii) the spatial distribution of the conductance changes recorded over the 
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O2[001] molecule for the 0-0 stretch and the antisymmetric 02-Ag stretch vibrations, (iii) 
the absence of signal for the center-of-mass and hindered rotation modes, and (iv) the lack 
of lETS signal for the molecule chemisorbed along the [110] direction. These results permit 
us to give an enhanced description of the JETS signals in terms of the molecular symme- 
tries of the adsorbed system. Leading to the conclusion that a multi-orbital description is 
mandatory and that the simulation is able to capture the complex increases and decreases 
of conductance depending on the STM tip location over the molecule. 

The good agreement of the simulations with the experiment validates the results of the 
electronic structure as computed with DFT+f/. In particular, DFT+f/ shows that out of the 
two possible molecular conformations on the surface, only the one along the [001] direction 
is paramagnetic. 

The inclusion of the U correction does not alter the values of the other DFT parameters 
such as the molecular geometry and vibrational structure. Regarding vibrations, in order to 
attain reliable numerical values of the molecular frequencies, we have used a Morse potential 
to fit the DFT potential energy surface because anharmonicities rapidly set in and good 
frequency values are difficult to obtain. 

Finally, Fano-like profiles are revealed for the lETS as a function of the bias voltage, how- 
ever the two main sources of experimental broadening (temperature and more importantly 
the lock-in modulation voltage) erases any information of the actual spectral line shapes. 
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TABLE I: Mode frequencies (meV) obtained by DFT+U (?7=3 eV) and diagonalization of the 
dynamical matrix computed with the same displacement AR for all atomic oxygen coordinates 
with O2 on Ag(llO) along the [001] direction. The modes are schematically depicted in Fig. [5j 

Mode AR = 0.025 A AR = 0.03 A Ai? = 0.04 A Experiment 
0-0 stretch 91.5 89.7 90.4 82.0^4 
antisym. 02-Ag 42.0 31.8 33.6 38.aP 
center-of-mass 35.2 22.5 25.1 3CP™ 
hindered rotation 34.5 37.8 34.5 



TABLE IL GGA and DFTH-?7 (?7=3 eV) quantities describing the physical properties of the two 
chemisorption states of O2 on Ag(llO). 










0-Ag 


02-surface 


magnetization 


state 


calculation 


(A) 


(A) 


(A) 




O2[001] 


GGA 


1.43 


2.35 


1.34 


0.0 




GGA + U = 3eY 


1.41 


2.37 


1.34 


0.54 


O2[110] 


GGA 


1.45 


2.36 


1.12 


0.0 




GGA + [/ = 3 eV 


1.48 


2.35 


1.12 


0.0 
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FIG. 1: Potential energy curve along the 0-0 stretch mode of O2 on Ag(llO) (squares) and a 
Morse potential fit (black line). The obtained frequency is 83.6 meV in excellent agreement with 
the experimental one (82.0 meV^. 
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FIG. 2: (Color online) Spin-up (positive) and spin-down (negative) projected density of states 
of O2[001] (left panels) and O2[li0] (right panels) onto the O2 molecular orbitals. (a) and (c) 
Semilocal DFT (GGA) results, (b)and (d) Semilocal DFT+t/ results (U =3 eV). Intra-atomic 
correlation leads to a change in the TVg resonance for the O2[001] only. The Gaussian broadening 
used in the projections is 0.25 eV. 
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FIG. 3: (Color online) Calculated DFT+t/ {U =3 eV) local DOS for O2 on Ag(llO) calculated 
with a sample bias V=200 mV. The zero height corresponds to a tip-surface distance of ^ 6 A. 
(a) and (b) Topographical image and profile (both in A) along the molecular axis for O2[001]. The 
protrusions attributed to the tt^ orbital are indicated by arrows, (c) and (d) Topographic image 
and profile (both in A) along the molecular axis for O2[li0]. The white grid lines in (a) and (c) 
show the surface structure. 
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FIG. 4: (Color online) Calculated DFT+t/ {U =3 eV) STM-IETS for O2 chemisorbed on Ag(llO) 
along the [001] direction. Spatial distribution of the conductance changes for: (a) the symmetric 
0-0 stretch mode {V=200 mV) and (b)the antisymmetric 02-Ag stretch mode (y=150 mV). (d) 
and (f) Cross sections of (c) and (e), respectively, along the molecular axis. Note the negative 
and positive conductance changes observed on the antisymmetric mode, in contrast to the overall 
negative values of the 0-0 stretch mode. The white grid lines show the surface structure. 
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FIG. 5: (Color online) Sketch and symmetries of the MOs and modes involved in the lETS images. 
The table summarizes the MOs that can couple with the vrj- to assure a nonzero matrix element 
('0m,k|'^|'0n,k) [scc Eqs. and Q]. The symmetric (S) or antisymmetric (A) character of the tv^ 
and the electron-vibration coupling v is written in black following this order for each vibrational 
mode and symmetry plane. In gray (red), the symmetry character of '0m,k necessary for a nonzero 
coupling. The MOs fulfilling those symmetry conditions are written inside the gray (red) squares. 
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FIG. 6: (Color online) Contributions to the conductance changes shown in Fig. |4]For O2[001]. (a) 
Inelastic and (b) elastic contributions to the 0-0 stretch mode, (c) Inelastic and (d) elastic to the 
antisymmetric 02-Ag stretch mode. 
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FIG. 7: (Color online) Contributions to the conductance changes for O2[001]. (a) Inelastic and 
(b) elastic contributions for the center-of-mass stretch mode, (c) Inelastic and (d) elastic for the 
hindered rotation mode. 
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V(V) 

FIG. 8: (a) Conductance of O2 on Ag (110) along the [001] direction in a model tunnel junction. 
The 0-0 stretch mode is allowed to interact with the tunneling electrons leading to a decrease of 
conductance over the vibrational threshold, (b) Derivative of the conductance with respect to the 
applied bias. Both graphs show the behavior with temperature. In (b) the full line spectra have 
been convoluted with a gaussian broadening representing the effect of a 7 mV rms modulation 
voltage, (a) and dashed lines in (b) are without modulation-voltage broadening. 
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